Fractal Properties and Small-scale Structure of Cosmic String Networks 
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We present results from a detailed numerical study of the small-scale and loop production prop- 
erties of cosmic string networks, based on the largest and highest resolution string simulations to 
date. We investigate the non-trivial fractal properties of cosmic strings, in particular, the fractal 
dimension and renormalised string mass per unit length, and we also study velocity correlations. We 
demonstrate important differences between string networks in flat (Minkowski) spacetime and the 
two very similar expanding cases. For high resolution matter era network simulations, we provide 
strong evidence that small-scale structure has converged to 'scaling' on all dynamical lengthscales, 
without the need for other radiative damping mechanisms. We also discuss preliminary evidence 
that the dominant loop production size is also approaching scaling. 



I. INTRODUCTION 

Topological defects are an unavoidable consequence of 
phase transitions in the early universe (for a review see 
ref. I'd. Cosmic strings, in particular, will form in a 
range of cosmological scenarios, including at the end of 
brane inflation 0^0, 3 . Since they are intrinsically non- 
linear objects, one is ultimately forced to resort to high 
resolution numerical simulations 0, if one wants to 
study their evolution and cosmological consequences in 
detail. 

In recent years, while some authors have exploited 
ever-increasing computing capabilities 0, 0, 0j many 
others have tried to simplify the problem. For example, a 
number of flat (Minkowski) spacetime string simulations 
have been performed (lOl, [ul El E El El , which it has 
been argued are a good approximation to the expand- 
ing case. (The idea is that the string correlation length 
is significantly smaller than the horizon size and so the 
network should not feel the expansion.) Another popular 
approach is to generate artificial networks [la. Il7j made 
of a (variable) number of segments with a given size and 
velocity, and aiming to mimic the evolution of a real net- 
work by enforcing (by hand) suitable variations of these 
parameters (e.g., as predicted by quantitative analytic 
string evolution models E1>ESI13)- Finally, doubts 
have also been raised as to whether the Nambu action 
accurately reproduces the dynamics of vortex-strings in 
the underlying field theory 0, (though we will not 
discuss this issue further here, see Q). The problem 
is that different approaches yield quantitatively differ- 
ent and even sometimes contradictory cosmological con- 
sequences. This raises the question of how reliable all 
these approximations really are. One of the goals of this 
letter is to clarify these issues. In particular, we argue 
that there are important differences between string net- 
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work evolution in flat spacetime and in an expanding 
universe. 

There is also considerable confusion in the literature 
on the issue of scaling. This is partially due to the fact 
that people define scaling in different ways. The weak- 
est definition ~ large-scale scaling - is simply to say that 
one wants the energy density to be a constant fraction 
of the total density in the universe. This is the attractor 
solution of the evolution of long-string networks, and is 
strongly confirmed by both analytical and numerical ev- 
idence. On the other hand, this by itself does not guar- 
antee that small scale features of the network are also 
scaling, such as the typical scale of small scale wiggles, 
that of loop production, or even the correlation length 
itself. In this letter we aim to clarify these different con- 
cepts. We also present strong evidence for scaling of the 
small-scale features of the string network (as measured 
by its multifractal dimension [231 123. l25l| and renormal- 
ized mass per unit length) and we discuss the scaling of 
the characteristic loop production size. 



II. AVERAGED AND MICROSCOPIC 
PROPERTIES 

We have performed many thousands of CPU hours of 
cosmic string simulations, using an upgraded version of 
the AUen-Shellard string code on the COSMOS super- 
computer. Some preliminary, large-scale results have al- 
ready been discussed in and a thorough description 
of the simulations will be presented elsewhere Here, 
however, we take advantage of the extremely high reso- 
lution of our runs to provide a detailed description of the 
averaged and small-scale properties of cosmic string net- 
works, as well as of the processes that determine loop 
production. This goes substantially further than the 
previous analyses of refs. [27L |2^ and it is interesting 
to compare with the more recent flat spacetime analysis 
of ref. 15j. Ultra-high resolution simulations were per- 
formed in the matter and radiation epochs and in flat 
(Minkowski) spacetime. The initial Vachaspati-Vilenkin 
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networks '2^ have resolutions between 75 and 256 points 
per correlation length (PPCL), and we subsequently en- 
force a constant resolution in physical coordinates. This 
resolution is a substantial improvement over currently 
published simulations. These networks are then evolved 
for a dynamic range (in conformal time) of order 3 (and 
up to 6), each taking several weeks of processor time. 

We note that not everybody agrees on when a network 
reaches the scaling regime. For example, numerically we 
find that the number of strings per horizon pt^/^Q be- 
comes a constant (a plausible definition of linear scal- 
ing) before the correlation length ^ becomes a constant 
fraction of physical time t (which is another possibil- 
ity). This, in turn, happens before other properties of 
the network reach scaling — indeed, there are a number of 
different timescales involved in the approach to scaling. 
We have developed our own statistics, based entirely on 
small-scale properties, to quantify how far a network is 
from scaling. We will describe them in detail elsewhere 
|26j. but for the purposes of this letter it is enough to 
think of the two definitions above as roughly measuring 
large-scale and small-scale scaling. Our high resolution 
and large dynamic range allow us to go well beyond the 
point where large-scale scaling is reached, and even to 
reach small-scale scaling (under certain circumstances). 

We must emphasise that we have chosen initial con- 
ditions for the radiation and matter era boxes for which 
we knew (by performing other smaller simulations be- 
forehand 's',^ and by using analytic modelling (2(1. l2ll|) 
that the approach to scaling would be fastest. Also with 
this aim, we have started the runs giving the strings an 
initial RMS velocity Vi ~ 0.6. For other choices of initial 
conditions, the approach to scaling would be consider- 
ably slower. By the end of each run, the error made in 
energy conservation was less than half of a percent. 

Consider first some averaged properties of the net- 
works, displayed in Fig. ^ This shows the evolution of 
the 'characteristic' lengthscale 

Poo=Mo/i^, (1) 

with ^0 being the 'bare' string mass per unit length), 
the correlation (or persistence) length ^ (to be defined 
below) , the RMS string velocity, the coherent velocity at 
the scale of the correlation length Vc, the dimensionless 
parameter poct^ /p-o (essentially a measure of the num- 
ber of long strings inside one horizon volume), and the 
'renormalised' dimensionless string mass per unit length 
defined as 

/Ltoff = p.p.0 (2) 

at the scale of the correlation length. Note that the fact 
that we have started the runs with a non-zero string ve- 
locity is responsible for giving the initial network an effec- 
tive mass per unit length different from unity. Note also 
that while L is simply just a measure of the string energy 
in a given volume, ^ is the actual correlation length, and 
it is measured directly from the simulation. It is obvious 
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FIG. 1: The time evolution of some characteristic large-scale 
properties of cosmic string networks, for 75 PPCL matter 
(solid lines), radiation (dashed) and flat-spacetime (dotted) 
runs. In all plots the horizontal axis represents the dynamic 
range in conformal time. All plotted quantities are defined in 
the text. 

these properties have strikingly different behaviours in 
the three different regimes. Interestingly, the only prop- 
erty which seems to be relatively independent of the cos- 
mological scenario is the coherent velocity at the scale of 
the correlation length. 

We have also performed a detailed analysis of the 
small-scale properties of these networks. This involves a 
significant amount of data analysis, the details of which 
will be provided elsewhere (26i |. Here, we will simply anal- 
yse a typical 'late-time' box of each of the three simula- 
tions above. (The analysis of the time evolution of these 
properties for different epochs in the approach to scaling 
can also be done.) In order to have a more meaningful 
comparison, we have chosen boxes that have evolved for a 
dynamic range (in conformal time) just above 3 (the lat- 
est time-step we obtained for the radiation era run). We 
have exhaustively analysed boxes for other nearby time- 
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FIG. 2: Some characteristic small-scale properties of cosmic 
string networks in the hnear scahng regime, for 75 PPCL mat- 
ter (sohd hnes), radiation (dashed) and flat-spacetime (dot- 
ted) runs. In all cases, the properties are measured after evo- 
lution by a dynamic range of about 3. In all plots the horizon- 
tal axis represents the logarithm of the physical lengthscale 
relative to the correlation length of the network. All plotted 
quantities are defined in the text. 



steps, and find that the results to be presented below 
remain unchanged, within errors, so we can confidently 
say that these properties are indeed representative of the 
epoch in question (with one exception that we will point 
out below). The results are presented in Fig |21 as a 
number of quantities measured on a range of physical 
scales relative to the correlation length ^. Each point on 
the graphs results from an average over (typically) a few 
thousand measurements performed on randomly selected 
points on the simulation box. The plots do not show 
the entire range of scales we probe — only those which we 
believe to be accurately sampled and free of systematic 
errors (e.g., box-size effects). 

The top two panels show the correlation function for 



the tangent vectors, 

corrx{s) = (x'(so) • x'(so + s)) (3) 

and its cumulative sum {Cx{s) = corrx{s')ds' , plotted 
relative to physical time t). Notice the differences be- 
tween the three cases for scales around and sightly below 
the correlation length. In particular, this confirms that 
the correlation length ^ (which is defined as the smallest 
scale at which the correlator vanishes, and is therefore 
roughly the height of the peak in the Cx/t curve) is quite 
different in the three cases, as we already saw in Fig. ^ 
The middle panels show the analogous functions for the 
velocity vectors, x. A first striking feature is that in 
the expanding universe cosmic string velocities are anti- 
correlated on scales between the correlation length and 
the horizon. However, such a feature is not present in fiat 
spacetime. This anti-correlation is the result of a 'mem- 
ory' of the network for recent intercommutings, and its 
absence in the fiat spacetime case highlights the fact that 
the loop production mechanism is different in the expand- 
ing and non-expanding cases. Indeed, such an effect was 
discussed in (see also 113). If one defines a 'velocity 
coherence length', this will be significantly smaller than 
^ itself. 

Finally, the bottom two panels show what we believe 
to be the most fundamental microscopic properties of the 
string networks. On the left we have plotted the general- 
ized fractal (also referred to as 'multi-fractal' or 'differen- 
tial fractal', depending on the context) dimension of the 
network |23l 1241 l25l | , again as a function of scale relative 
to the correlation length. On the right we have the 'renor- 
malized' (dimensionless) string mass per unit length, ^. 
The fractal dimension is obviously unity on small scales 
(strings are roughly straight lines) and two on very large 
scales (strings are roughly random walks). The inter- 
esting question, however, is what happens at intermedi- 
ate scales. In particular, one should expect a range of 
scales where strings should behave as self-avoiding ran- 
dom walks, and these have a fractal dimension d^ — 5/S 
in three spatial dimensions. 

We can immediately see that the radiation and mat- 
ter era boxes have strikingly similar fractal profiles (re- 
call that we have rescaled by the respective correlation 
lengths, so the physical scales involved will be different). 
The differences result in a slightly higher mass per unit 
length on a given scale (as one can confirm on the right). 
The flat spacetime profile, however, is qualitatively differ- 
ent, and shows a rather high 'step' around the scale of the 
correlation length. We have evidence that this step in the 
fractal dimension tends to diminish and move very slowly 
to smaller scales as the simulation progresses, indicating 
that fiat space simulations have not really achieved small- 
scale structure scaling. Nevertheless, we have performed 
moderately high resolution simulations with a dynamic 
range of 10 and the 'step' is still clearly visible, so we 
believe it is to some extent a persistent feature. This dif- 
ference in the fractal properties must refiect the differing 
efficiencies of loop production in fiat space and the ex- 
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panding universe. While the integrated loop production 
efficiency is much greater in flat spacetime c = 0.57, as 
opposed to the expanding c = O.23]370j it appears to be 
relatively less effective around the correlation scale with 
energy 'stuck' on fairly large scales. This intuitive pic- 
ture can be confirmed by noting that the renormalised 
mass per unit length /i is smaller than in the expanding 
case on small scales, but is larger on large scales. 

Relevant to this issue, we note that a different (func- 
tional forms) approach to flat spacetime simulations was 
recently presented by |15| . Simple comparisons indicate 
that their results appear to be in good agreement with 
our flat spacetime ones, though note that an allowance 
has to be made for the fact that the two codes use dif- 
ferent definitions when calculating the correlation length 
^. Another difference is that our characterization of the 
statistical properties of the network is based on fractal 
properties, whereas theirs is based on Fourier space tech- 
niques. In fact, there are very interesting relations be- 
tween the two approaches, that will be explored elsewhere 

Finally, we point out that in all cases the fractal di- 
mension of the network at the scale of the correlation 
length is well below two: typical values are 1.2 in the ex- 
panding case and 1.4 in the flat case. (Note that this 
is to be expected if one interprets ^ as a persistence 
length.) So the intuitive picture that a string network 
looks Brownian at the scale of the correlation length is 
clearly incorrect. Indeed, it's not even strictly true at the 
scale of the horizon — here the network looks more like a 
self-avoiding random walk (an obvious consequence of in- 
tercommutings) . The Brownian picture is only valid on 
significantly larger scales. This is relevant in a number of 
contexts, in particular for the 'toy model' approach dis- 
cussed in the introduction, which implicitly models the 
network under the former assumption. 



III. WHEN ARE WE SCALING? 

The key question is now whether one can see scaling 
of the small-scale features of the cosmic string network. 
Detailed analysis indicates that quantities like the mul- 
tifractal dimension dm, renormalized mass fj, or coherent 
velocity as a function of scale evolve significantly at the 
beginning of the simulations, and even continue to do 
so after the large-scale properties of the network reach 
scaling (say, as measured by L/t = const, so that there 
in a constant number of strings — about 40 — per Hubble 
volume). On the other hand, if one has high enough res- 
olution (to probe lengthscales well below the correlation 
length) as well as dynamical range, one will see the evo- 
lution of dm{() and /i(£) stop, hence achieving scaling. 
The first evidence for this is shown in the left panels of 
Fig. where these quantities are plotted for the late-time 
evolution of five different matter era runs. Note that the 
initial conditions in each are very different. Not only do 
the resolutions vary between 75 and 256 PPCL, but the 



128 and 256 PPCL boxes are two very different simula- 
tions, one with Vi — and the other with vt = 0.6 (and 
with different initial densities). Despite this, one sees 
that the dm{£) and /i(^) profiles are virtually identical, 
within errors. We have similar evidence of scaling for 
the radiation era, although the convergence is somewhat 
slower — this will be discussed in more detail elsewhere 

These matter era results, in particular, now offer 
unambiguous evidence that small-scale structure has 
achieved scaling on all lengthscales accessible to these 
high resolution numerical simulations. As indicated by 
and conjectured from an earlier study 27], it now seems 
clear that the dynamical processes of reconnection and 
loop production are sufficient to govern the build-up of 
small-scale structure on cosmic string networks, without 
requiring additional damping mechanisms such as gravi- 
tational radiation. Of course, simulations of even higher 
resolution and longer dynamical range will inevitably ac- 
quire small-scale structure on even smaller scales, but it 
seems reasonable to suppose, given the convergence of the 
renormalised tension, that the integrated effect of these 
additional modes will be relatively small. We note that 
similar evidence for the scaling of small-scale structure 
on particular lengthscales has also been presented from 
flat spacetime simulations |0, though this is somewhat 
weaker because, like our own flat space results, a wide 
range of other lengthscales continue to evolve slowly in 
this case. 

A related but different issue is that of loop produc- 
tion. On the right-hand panels of Fig. |21 we quantify 
the net loop production function of each of these net- 
works for the timesteps corresponding to the left-hand 
plots. We show both the overall net production func- 
tion (top panel, summing all loop creating or destroy- 
ing events) and the production from long strings only 
(bottom panel). The dominant overall loop production 
mechanism is that of the self-intersection of larger loops, 
so note that direct loop production by long strings typi- 
cally happens on larger lengthscales (above the resolution 
of the simulations). It is important to note that these are 
the time-dependent loop production functions, that is, 
the loops produced or destroyed on each lengthscale in a 
very narrow time span centered around the corresponding 
timestep. This is a much finer diagnostic of loop genera- 
tion mechanisms than the overall string loop density, and 
it also indicates the typical lengthscales of the long-lived 
loops which are relevant, for example, for gravitational 
wave background estimates. Unfortunately, Fig. |2| indi- 
cates that the dominant loop production scale remains 
at least weakly dependent on the simulation resolution - 
there is not the same strong agreement as between the 
fractal profiles on the left-hand side. This is even clearer 
from the equally-spaced time evolution of the loop pro- 
duction function shown in Fig. 0] for the 75 PPCL mat- 
ter era. This scale starts out being about the size of the 
correlation length, but becomes progressively smaller as 
small-scale structure builds up on the strings. The evo- 
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FIG. 3; The multi-fractal dimension (top left), renormalized 
mass per unit length (bottom left), overall net loop produc- 
tion function (top right) and net loop production function 
from long strings only (bottom right) for a series of 5 different 
matter era runs. The resolutions (in PPCL) and dynamical 
ranges (in conformal time) are indicated in the legend. 



lution of the peak of the loop distribution, however, is 
clearly beginning to slow down at late times indicating 
that it is rising above the minimum simulation resolution 
and will approach scaling. 

It is clear from Fig. |3| that even though all the net- 
works are statistically almost identical on a wide range 
of lengthscales, the dominant loop production scales are 
much more sensitive to variations on very small scales. 
Indeed we find that even when small-scale structure has 
reached scaling (as measured by dm{£), /^(^) or other di- 
agnostics to be discussed in the typical loop pro- 
duction scale can still evolve significantly towards smaller 
scales. However, we do have some evidence that this mi- 
gration is significantly slowing down, particularly in the 
matter era, so that the loop production scale will itself 
reach scaling. 
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FIG. 4: The overall net loop production function for differ- 
ent times in the evolution of the 75 PPCL matter era run. 
The dynamical range of each plot, equally spaced in terms of 
horizon times, are indicated in the legend. 



IV. DISCUSSION AND CONCLUSIONS 

We have presented results from a detailed analysis of 
the small-scale and loop production properties of cosmic 
string networks, in both flat spacetime and the expanding 
universe. We have shown that, while the differences be- 
tween the radiation and matter epochs are largely 'quan- 
titative', the evolution in flat spacetime is fundamentally 
different from the expanding case. This indicates that flat 
spacetime simulations, while qualitatively useful, will not 
provide an accurate approximation for cosmic strings in 
an expanding universe. 

The two distinguishing characteristics of string evolu- 
tion in Minkowski spacetime are the absence of velocity 
anti-correlations on scales around the correlation length, 
and the apparent existence of a 'preferred' scale (roughly 
around the correlation length ^) from which energy does 
not move to smaller scales. We believe that these can 
have a substantial influence, for example when calculat- 
ing unequal time correlators for the network, which we 
are currently investigating in more detail. Another result 
of our detailed analysis is that cosmic string networks are 
only random walks on length scales well outside the hori- 
zon, indicating that the semi-analytic models for strings 
need to be improved. 

Finally, from high resolution matter era simulations we 
have presented comprehensive evidence over all length- 
scales that small-scale structure converges to a scaling 
solution, that is, through reconnections alone without the 
need for additional damping mechanisms such as gravi- 
tational radiation. Evidence for similar convergence has 
also been found from radiation era simulations. While 
loop production is clearly more sensitive to the presence 
of structure on the smallest scales probed by the sim- 
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ulations, there are good indications that this will also 
achieve scaling. We will present a more detailed and 
quantified analysis of these results elsewhere . 
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Note added: While our paper was being completed, 
ref. [son] appeared on the archive. This paper presents 
some evidence for the scaling of the overall loop distri- 
bution on intermediate length scales below the correla- 
tion length but still near it, roughly loop lengths ^/lO 
- ^. A comparison with the net loop production func- 
tion (Fig. indicates that a population of short-lived 
relatively large loops is being characterised. These loops 
will soon cascade through self-intersection down to much 
smaller sizes, presumably determined at present by their 
simulation resolution. 
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